On the description of two-particle transfer in superfluid systems 
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Exact results of pair transfer probabilities for the Richardson model with equidistant or random 
level spacing are presented. The results are then compared either to particle-particle random phase 
approximation (ppRPA) in the normal phase or quasi-particle random phase approximation (QRPA) 
in the superfluid phase. We show that both ppRPA and QRPA are globally well reproducing the 
exact case although some differences are seen in the superfluid case. In particular the QRPA 
overestimates the pair transfer probabilities to excited states in the vicinity of the normal-superfluid 
phase transition, which might explain the difficult in detecting collective pairing phenomena as for 
example the Giant Pairing Vibration. The shortcoming of QRPA can be traced back to the breaking 
of particle number that is used to incorporate pairing. A method, based on direct diagonalization 
of the Hamiltonian in the space of two quasi-particle projected onto good particle number is shown 
to improve the description of pair transfer probabilities in superfluid systems. 

PACS numbers: 25.60. Je,25.40.Hs,24.30.Cz,21. 10.Re,21. 60. Jz 
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I. INTRODUCTION 

The importance of pairing correlations in nuclear sys- 
tems has been established in several aspects, binding 
energies of nuclei, odd-even effects, superfluid phenom- 
ena and pair transfer mechanisms, to mention just a few. 
However, despite the fact that pairing is anticipated to 
play a significant role in the pair transfer process, the 
existence of collective effects, like Giant Pair Vibration 

(GPV) pHl 

leading to an increase in the pair transfer 
from a superfluid nuclei still challenges the experimen- 
tal nuclear physics @- On the theoretical side, mean- 
field methods based on Hartree-Fock-Bogolyubov (HFB) 
sometimes augmented by Quasi-particle Random Phase 
Approximation (QRPA) have been used to predict pair 
transfer probabilities either from ground state to ground 
state or from ground state to excited states |cHL2|. A 
common conclusion of most of these studies is the sen- 
sitivity of one- or two-nucleon transfer process to the 
internal topology of pairing in nuclei. In the present 
work, exact results of pair transfer probability are ob- 
tained for the Richardson model [13J consisting in set of 
single-particle levels interacting through a pure pairing 
interaction. This model can be seen as the valence space 
of the last occupied level in nuclei where nucleons can be 
either added (pick-up reactions) or removed (stripping 
reactions). The possibility to perform exact calculations 
open new perspectives to understand the pair transfer 
process and provides a benchmark for approximate treat- 
ments. In the following, we first discuss the pair transfer 
mechanism from a general point of view and estimate 
pair transfer probabilities in the pairing model. 

We are interested here in a process where two particles 
are either added or removed in a system that is initially in 
its ground state formed of N particles. In the following, 
we denote by \v,A) respectively the eigenstates of the 
systems with A particles associated to the set of energies 
and, by convention, v = is taken for ground state. 
During its evolution, the system wave-function can be 



decomposed as [3] 
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where the first line describes the possibility that the sys- 
tem remains in its ground state or in one of its excited 
states without changing its particle number. Second 
(resp. third) line contains the information on the re- 
moval and/or addition process. The explicit form of the 
coefficients depends on the physical process under in- 
terest, like the stripping or pick-up reactions in nuclear 
physics. On the theoretical side, information of these 
processes can be obtained by studying the small ampli- 
tude response of the system to an external field T that 
changes the particle number by two units. Then, infor- 
mation on the transfer reduces to the knowledge of the 
response function, given by: 

S(E) = ^ \( N + 2 > v\f\N,0)\ 2 S (E - AE^ +2 ) 

+ J2\( N - 2 ^\T\N,0)\ 2 S{E-AE^- 2 ) (1) 

= S Add (E)+S Rcm (E) 

where AE^ ±2 = E^ ±2 - Eg. In the following, S Add (E) 
(resp. S Rem (E)) will be referred to addition (removal) 
strength function. A common choice of T @ to excite 
pairing modes is 



(2) 



where ct\at corresponds to creation operators of a pair 
of time-reversed single-particle states. In the present ar- 
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ticle, we are interested in a physical process where two- 
particles are either added or removed. In that case, it is 
more suited to consider directly the non-Hermitian ad- 
dition (resp. removal) transition operator, denoted by 
<f Add ( regp yRcm) de fined through 



f Add = J2 T u a H T Rem = ^2 T i 



gOjOi (3) 



where the particle-number operator Ni and pair cre- 
ation/annihilation operators PJ, Pi are given by 



Ni = a\a,i + a]ai , 



H = 44 



P = {P})^ (7) 



respectively associated to S Add (E) and S Rem (E). Prom 
the expression of the strength, we see that the under- 
standing of the two-particle transfer passes through a 
good knowledge of the spectroscopy of initial and final 
states as well as of the capacity to provide the addition 
or removal probabilities defined here as: 



P Add = |(A + 2,i/|T|A,0)| 2 , 

pRem = |(AT_2,Z/|f|7V,0)| 2 . 



(4) 
(5) 



If the many-body problem can be solved exactly, such 
quantities as well as the exact eigenvalues of the Hamil- 
tonian can be used to have a precise estimate of the 
strength function (JTJ. In most realistic situations, such 
treatment is impossible and educated guess should be 
used. In the present work, we are interested in an initial 
system where pairing correlation plays a role. This case 
has been first considered in refs. [if, [H, OH leading to 
the concept of pairing vibration, i.e. a coherent excita- 
tion of pairs of particles that is expected to show up in 
the enhancement of pair transfer probabilities. The stan- 
dard way to incorporate pairing correlation is to use the 
BCS or HFB approach as a starting point [ijlBll- This 
technique is indeed standardly used in nuclear physics ei- 
ther to estimate the transfer from ground state to ground 
state d, [H, [lj| or from ground state to excited states. 
In the latter case, the response is obtained using QRPA 
pi IH [pO] or its time-dependent version H. For a com- 
prehensive introduction, please refer to [201 ] . 

With the increase of computational powers, it is possi- 
ble nowadays to study exactly pair transfer in schematic 
model that approaches realistic situations and to quan- 
tify the predictive power of mean- field based approaches. 
In the present work, we first present exact results of 
pair transfer probabilities for the Richardson model with 
equidistant or random level spacing. The exact results 
are then used to benchmark standard approaches, namely 
ppRPA and QRPA. 



II. EXACT DESCRIPTION IN THE 
RICHARDSON MODEL 

A system of SI doubly-degenerated single-particle lev- 
els interacting via a pairing force with parameter G is 
considered. The Hamiltonian is given by [l3T | 



1 

0.5 


1 

■1 o 

i 1 

^0.5 
P 

1/5 0.4 

0.2 


0.4 

0.2 





_ (a) 








I'll 


(b) 


Exact 

ppRPA- 

\ ~ QRPA 


I'll 






i'ii 


- (d) /\ 


( ' /v N 


I'll 


(e) 


hJ v/V\ 







E (MeV) 



10 12 14 16 



H 



(6) 



FIG. 1: (color online) Addition strength function as a func- 
tion of the excitation energy E for different pairing strength 
G obtained using exact eigenvalues and eigenstates in Eq. (fl]) 
(red solid line) for a system of iV = 8 particles to N = 10 in 
the case of — 10 equally spaced levels. From top to bot- 
tom: (a) G/Ae = 0.1 , (b) 0.3, (c) 0.5, (d) 0.7 and (e) 0.9 
are shown. The results of the ppRPA (blue dashed line) and 
QRPA (black dot-dashed line) are also presented. Note that 
in panel (c), the result of ppRPA obtained at the collapse 
point G = 0.48Ae is also reported. The QRPA are shown 
only above the BCS threshold G = 0.33Ae. 

For not too large model space fi, exact solutions can 
be obtained using standard diagonalization techniques in 
subspace of given seniority [21( . As a test case here, a sys- 
tem of = 8 particles with SI = 10 doubly degenerated 
levels is considered here. Assuming a simple transition 
operator @ with = 1 for all pairs, illustrations of 
addition strength function obtained for the Richardson 
model are shown in Fig. [TJ in full (red) line, for equidis- 
tant level spacing, i.e. ti = iAs (i = 1, SI), for different 
values of the pairing interaction and Ae = 1 MeV. In the 
following the excitation energies will be calculated with 
respect to the ground state of the system with N = 10 
particles consistently in all the theories. The presented 
result is exact in the sense that the exact eigenvalues 
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and eigenstates of the system with N = 8 and N = 10 
have been used to compute Eq. ((T|). Note that in the 
present work, we are mainly interested in the transition 
from ground state to excited states of the N + 2 nu- 
cleus and the contribution of the ground state to ground 
state transfer has been omitted in the figure. Moreover, 
in order to make simpler the comparison between differ- 
ent results, we have folded the discrete spectra with a 
Lorentzian function with a width of 1 MeV. For com- 
pleteness, we also show in Fig. [5] similar results obtained 
with randomly spaced system whose single-particle en- 
ergy are 0.551, 2.176, 4.033, 5.142, 6.444, 7.029, 7.827, 
8.343, 9.226 and 9.571 in MeV units. This situation could 
be considered closer to realistic cases. 



III. PPRPA VS QRPA APPROACHES TO PAIR 
TRANSFER 



As already mentioned, in most cases, exact evaluation 
of the pair transfer probabilities cannot be performed and 
approximations for the many-body states are necessary. 
The most common strategy used in nuclear physics is 
to first apply the HFB or BCS theory and minimize the 
energy in the Hilbert space of quasi-particle vacuum im- 
posing a mean particle number equal to N. This leads to 
an approximation for |0,iV). Using standard notations, 
the quasi-particle vacuum is given by 



|0,AT) ~ |0,QP) =IX [' 



a] at 



h) 



(8) 
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FIG. 2: (color online) Same as figure [T] obtained for randomly 
spaced single-particle energies. See the text for more details. 
The QRPA are shown only above the BCS threshold G = 0.44 
MeV. 



In the exact calculations, we observe a shift of the ex- 
citation spectrum toward higher energies as the pairing 
strength G increases while the pair transfer probabilities 
of the excited states get smaller. On the contrary, as it 
will be shown below, the pair transfer probability from 
ground state to ground state increases as G increases. 



where |— ) is the quasi-particle vacuum and a$ is the 
quasi-particle annihilation operator defined through the 
simple Bogolyubov transformation: 



a, = Uidi - ViaA. 
oca = UiOj+Via]. 



(9) 
(10) 



In the HFB approach, the above transformation auto- 
matically implies that the single-particle basis identifies 
with the canonical basis. 

In the present model, the Hamiltonian ^ is already 
written in the canonical basis for the HFB theory. This 
theory is particularly suitable to provide estimates for 
ground state (GS) to GS transfer probabilities [H, [Hj]. 
Since we want to compare with exact results, at the 
BCS/HFB and QRPA level the contribution to the 
particle-hole channel of the pairing interaction is taken 
into account. The result of the HFB theory for this prob- 
ability is shown in figure [3] and compared to the exact 
solution. Note that below the pairing threshold (denoted 
by G CI and equal to 0.33 MeV and 0.44 MeV for the 8 
particles system in the equally spaced and random spaced 
case, respectively), the HFB reduces to HF and probabil- 
ity is 1. As illustrated in this figure, while the HF theory 
(not shown) would have failed to reproduce the exact 
probabilities, the HFB framework gives estimations that 
are already rather close to the exact ones in the supcrfluid 
regime. 

Due to the absence of residual coupling between quasi- 
particle excitations, it is known that HFB alone can- 
not properly describe excited states. Then, linear re- 
sponse theory including possible particle-particle (pp), 
hole- hole (hh), or particle-hole (ph) excitations is applied 
to describe excited states, and then transfer probabilities, 
within the QRPA approach. In QRPA, the excited states, 
denoted by \v), are obtained by considering coherent su- 
perposition of 2 quasi-particle (2QP) excitations. This 
leads to 



W) = Ql\o), 



(ii) 



where Ql are QRPA phonons written as 



Ql=Y J {X u j a\a\-Y»cqa i l Q u = (Ql) f , (12) 
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while |0) is the phonon vacuum, defined through the con- 
ditions Q„|0) = 0. In practice, the components X v and 
Y v as well as the energies of the excited states ui v are 
deduced by solving the QRPA eigenvalue problem [22j . 
Since these techniques are rather standard [17j . we only 
recall here the expressions of the pair transfer probability. 
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FIG. 3: (color online) Illustration of the application of mean- 
field theory to provide estimate of ground state to ground 
state addition pair probability. The exact result (red circles), 
BCS (black triangles) and BCS projected on good particle 
number (blue open squares) are shown here for the equidistant 
single-particle level case. The result of the P-QTDA 1 - 08 - 1 is 
also presented with green filled diamond. 

In QRPA, the addition transition probability is given 
by 

p Add = | (Q|f Add^p = | J2(VM V) ~ UfY^)\\ (13) 

i 

It is well known that, the first solution of the QRPA 
equations corresponds to the spurious mode, which is 
then not considered in the evaluation of the strength 
function. 

In the weak coupling limit, below a certain threshold 
value of G denoted by G cr , the minimization of the en- 
ergy in HFB identifies to the Hartree-Fock approach with 
no pairing. Then the mean-field vacuum is a pure Slater 
determinant where the lowest hole states are occupied. In 
this case, labeling by h the hole state (V/, = 1, Uh = 0) 
and p (V p = 1, Up = 1) the particle states associated 
to this vacuum, the excited states are described by us- 
ing the particle-particle RPA (ppRPA) where the phonon 
creation operators (TT21 can be written as 

qI = £*p44+E y *" a /i4> (14) 

p h 

while the addition probability simply writes: 

= iE^-E^i 2 - as) 

p h 

Note that in the ppRPA case, contrary to the QRPA, 
the (7(1) symmetry associated to particle number con- 
servation is not broken. Explicit forms of the ppRPA 



and QRPA equations for the model considered here can 
be found in Refs. djEU 

In Figs [1] and [2j the ppRPA (dashed line) and QRPA 
(dot-dashed line) are compared to the exact results. 
Above a given threshold G^f , ppRPA collapses and 
leads to imaginary energies making not possible a direct 
comparison with the exact and QRPA results. However, 
in Fig. Q] we show in the panel (c) corresponding to a 
pairing strength G — 0.5MeV the ppRPA results ob- 
tained at the collapse point, i.e. G p r PA = 0.48, in order 
to show that its description is still reasonable even in the 
superfluid phase. 

From these comparisons, the following conclusions can 
be drawn: (i) The ppRPA does reproduce perfectly the 
exact results (energies and probabilities) in the normal 
phase, (ii) The QRPA provides a global reproduction of 
the pair transfer probabilities in the superfluid phase. In 
particular, the threshold in energy that is directly re- 
lated to pairing correlation is properly accounted for. 
It is worth to mention, that such a threshold can only 
be described in a mean-field theory by breaking particle 
number symmetry. Finally, it is also clearly seen that 
some differences exist between the exact and QRPA. in 
general QRPA leads to peaks in the strength that are 
at slightly higher energies compared to the exact solu- 
tion while probabilities are slightly overestimated. These 
differences are even stronger in the random space case 
and can stem from different origins. First, while part of 
the four quasi-particles (4QP) are accounted for in the 
QRPA ground state correlations, the complete inclusion 
of 4QP excitations is known to modify excited state en- 
ergy spectrum (l7l |26|. Second, a systematic error exists 
due to the breaking of particle number. Panel (c) of Fig. 
[1] illustrated that when ppRPA is applicable in the super- 
fluid phase G < 0.48 MeV, it gives a better agreement 
with the exact result compared to QRPA. Since ppRPA 
is a particle number conserving theory, this is a first indi- 
cation that the breaking of (7(1) symmetry might pollute 
the QRPA predictions. 



A. Role of particle number in the estimation of 
pair transfer probabilities 

At this stage, it is most likely to conjecture that the 
failure of QRPA to reproduce two-particle transfer pro- 
cesses stems from the mixing of systems with different 
particle numbers. The QRPA approach implicitly as- 
sumes that the states \v) are relatively good approxima- 
tions for the eigenstates of systems with N + 2 (or N — 2) 
particles. As an illustration of the correctness of this as- 
sumption, the mean number of particles N v = (y\N\v) is 
displayed in the bottom panel of Fig. @] as a function of 
the excitation energy w v . N v has been estimated using 
the quasi-boson approximation leading to: 

N U = J2 2(0? - V?)(X? + Y? 2 ) + (AO, (16) 
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where (TV) is the number of particle in the QP vacuum. 

In the same figure, the mean-particle number of the 
two quasi-particle (2QP) excited states \k) as a function 
of (efe — A) is also shown, where A is the Fermi energy. 
The 2QP states are defined through: 



|*) = alaU0,QP) 



(17) 



where the ground state has N particles in average. The 
mean particle number in \k) is given by: 



N k = (k\N\k) = (N)+2U 2 k -2Vl 



(18) 



This expression as well as the illustration in Figure 0] 
clearly shows that the 2QP states will be close to a state 
with TV + 2 (resp. N — 2) particles only if Uk — > 1, i.e. 
well above the Fermi energy A (resp. Uk — > 0, i.e. well 
below the Fermi energy), but will be a bad approxima- 
tion if the 2QP state involves single-particle state in the 
vicinity of A. Consequently, QRPA states will also suffer 
from the same problems if the state is constructed from 
2QP states that are close to the Fermi energy. While the 
QRPA results are in a reasonable agreement with the ex- 
act case, the effect of particle number conservation on the 
pair transfer is largely uncontrolled within QRPA. This 
is anticipated to be especially crucial in exotic nuclei as 
the level spacing is reduced closed to the drip line. 
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IV. IMPROVED TREATMENT OF PAIR 
TRANSFER IN SUPERFLUID SYSTEMS 

The effect of particle number conservation on pair 
transfer from ground state to ground state has been al- 
ready studied in Ref . [l2j . It has been empirically found 
that the breaking of U(l) symmetry has a rather small 
impact on the estimated probabilities. As a further il- 
lustration, we show in Figure [3] estimations of transfer 
probabilities with equation (j4} using the ground state 
quasi-particle vacua projected either on N or on N + 2 
particle numbers (see Ref. [l2j for technical details). As 
seen in Fig. [3l the BCS approach reduces to HF below 
the threshold and is not able to reproduce the trans- 
fer probability at low G. The projection after variation 
(open square) obviously does not cure this problem but 
considerably improves the probability above the thresh- 
old especially in the strong coupling regime. Note that 
As a conclusion, the HFB and/or BCS approach aug- 
mented by an a posteriori projection is already very good 
to describe GS to GS transition in the superfluid regime. 
Therefore, in the following we will concentrate the discus- 
sion on excited state where it is necessary to go beyond 
HFB. 

From previous discussion, we are facing the following 
dilemma: to describe the physics of pairing and in partic- 
ular the gap in energy between the ground state and the 
first excited state in a superfluid system, it is necessary 
to break the U(l) symmetry. On the other hand, this 
symmetry breaking seems to be at the origin of some dis- 
crepancies between QRPA and exact pair transfer prob- 
abilities. 

The most direct extension of the projection technique 
to estimate transition from ground state to ground state 
would be to directly estimate the transition from the 
BCS/HFB ground state projected onto N particle to the 
QRPA eigenstates projected onto N + 2 particles. This 
has however two disadvantages (i) in practice we found 
that the pair transfer probabilities that are much smaller 
than the exact ones (ii) QRPA states are not anymore 
orthogonal after projection leading to some difficulties to 
interpret the probabilities themselves. 

Alternatively, one can try to develop a RPA like ap- 
proach directly in the space of projected 2QP states. Fol- 
lowing the Tamm-Dancoff approximation spirit, a set of 
states j<E>fc) defined through: 



= P N+2 alal\Q,QP). 



(19) 



are introduced, where Pn+2 is the projector on N + 2 
particles [llj]. Then, excited states of the system with 
N + 2 particles decompose as 



FIG. 4: (color online) Bottom: mean number of particles N v 
of the QRPA states as a function of their excitation energy 
Twj v . Top: mean number of particles Nk of 2QP states as a 
function of (et — A), where A is the Fermi energy. 



\u,N+2)=J2Xk\*k) 



(20) 



This strategy has been analyzed in Ref. [27J as well as its 
RPA generalization following Refs. [28[ (see also [29l|). 
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A proper description would require a full Projected 
QRPA calculation whose practical implementation would 
be rather cumbersome, especially in realistic calculations. 
A simpler approach is to introduce what could be consid- 
ered as a projected version of a two quasiparticle Tamm- 
Dancoff approximation. Again, it should be mentioned 
that contrary to standard TDA, states \$k) are neither 
normalized nor orthogonal with each others. Therefore, 
a special attention has to be paid when formulating the 
approach. In practice, this implies to diagonalize the 
overlap matrix Oki = (&k\&i) prior to write the TDA 
eigenvalue problem. A practical method, has been pro- 
posed in ref. [27j to obtain the TDA equation in the 
projected space. 



Following 
written in terms of new states 



271 ] . the excited states are 



\u,N + 2) 



k 



States ($fc)) arc defined though: 

|(*fc)> = |$ fe )-{0 N +2|*fc), 



(21) 



(22) 



where |0n+2) correspond here to the approximated 
ground state with N + 2 particles. While in the original 
article, this ground state was anticipated to be obtained 
with a Variation After Projection (VAP) procedure, be- 
low it is simply taken as |0n+2) — P/v+2|0, QP). Then 
the TDA eigen-equation is solved in the space of |($fc)) 
states. In the following, this approach is referred to Pro- 
jected two Quasi-Particle TDA (P-QTDA). 

Consistently with the present approach, addition pair 
transfer probabilities are computed using the expecta- 
tion value of the transition operator between the quasi- 
particle vacuum projected on particle N and the excited 
states (ED: 



P 



Add 



Y,X%(Q,QP\P N fP N+2 \($ k )) 



(23) 



The projection onto different particle numbers in the ex- 
pression does not induce extra difficulty using the fact 
that 



(0,QP\P N TP N+2 \$ k ) 



(0,QP|TPtf +2 |$fc) 
(0,QP|JVf|$ fc ). 



Therefore expectation values entering in expression (|23j) 
can be performed using standard projections techniques. 
In practice, the projection on particle number is made 
by discretizing the gauge angle integration using the 
Fomenko approach with 199 points [lo, EH- Useful ex- 
pression to express the overlaps as well as the Hamilto- 
nian expectation value in a projected basis can be found 
in the appendix of Ref. [26] . 

Illustration of the method proposed in ref. [13] applied 
to pair transfer are shown in Fig. [5] and Fig. [6] respec- 
tively for the equidistant and random single-particle level 




E (MeV) 



FIG. 5: (color online) (color online) Same as figure [T] for 
(a) G = 0.5 MeV, (b) G = 0.7 MeV and (c) G = 0.9 MeV. 
The exact result (red solid line) is compared to the P-QTDA 
(black dashed line) and P-QTDA (GS) (blue dot-dashed line). 



spacing. These figures clearly demonstrate that the pro- 
jected TDA approach provides a much better description 
of the energy peak positions of the excited states in the 
system with N + 2 particles. However, probabilities of 
transfer are underestimated especially as G decreases. 

Further improvements can a priori be made by includ- 
ing more correlations in the ground state |0n+2)- Here, 
we simply used a Projection After Variation (PAV) ap- 
proach to approximate this state. A better treatment 
would be to perform a VAP as originally proposed in 
ref. Further correlations, like correlations induced 

by the coupling to 4QP states might be included using 
Projected QRPA instead of Projected TDA approach. 
The Projected QRPA approach has also been formulated 
previously (section III-B of Ref. [27]). In that case, not 
only the projected 2QP states but also the projected 4QP 
states should be explicitly introduced. 

The use of VAP and/or the introduction of 4QP, al- 
though possible in the present model (2(1 [HI, H3j|, will 
considerably increase the complexity of the approach in 
realistic situations (see for instance [3414361 ] for the VAP). 
Below, we propose a simpler method inspired from P- 
QTDA and able to grasp part of the ground state corre- 
lations without increasing the numerical complexity. 

Contrary to standard linear response theory based on 
HF (resp. HFB), the ground state projected onto good 
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FIG. 6: (color online) Same as figure[5]for the case of random 
level spacing. 



particle number is not orthogonal to the particle-hole 
(resp. 2QP) excited states. Similarly, the projected 2QP 
states are themselves not orthogonal to the projected 
4QP or higher order excitations. At first glance, this 
might be seen as a disadvantage compared to RPA or 
QRPA since additional orthonormalization is required. 
On the contrary, one might take advantage of the pres- 
ence of higher order components to improve the descrip- 
tion of the GS itself. 

The aim of the P-QTDA approach was to describe exci- 
tations with respect to the projected ground state. Then, 
the latter state has been naturally removed before solv- 
ing the eigenvalues equation (Eq. (|22p). Let us assume 
that we restart from expression (|2"0|) where the projected 
mean-field is also included in the sum ( with the conven- 
tion that |<J>o) = P N+2 \0,QP)). Then, coefficients X v k 
can be obtained by diagonalizing the hamiltonian in the 
{$fc} (with proper treatment of the non-orthonormality 
of the states). This direct approach, called below P- 
QTDA( GS ) , not only provides a way to get excited states 
but might also improve the description of the ground 
state itself. This is indeed what we observed empiri- 
cally. For instance, for G = 0.5 MeV and equidistant 
level spacing, the new ground state has an energy 300 keV 
lower than the energy of the original projected mean-field 
ground state. The difference reduces has G increase. At 
G = 0.9 MeV the reduction is only 30 keV. For complete- 
ness, the probability to transfer from GS to GS within 



the P-QTDA (GS) is also shown in Figure [3] This prob- 
ability is slightly improved compared to the BCS case 
below the BCS threshold while it follows the PAV results 
above. 

It turns out, that this approach improves the descrip- 
tion of pair transfer from GS to excited states. Results 
of the P-QTDA (GS) method are presented in Fig. (0 
and ([5]) by dashed line. A clear improvement is observed 
especially at low excitations and low G. The remaining 
difference with the exact solution is acceptable in view of 
the simplified approach presented here. 




0.4 0.6 
G (MeV) 



FIG. 7: (color online) Ratios of probabilities estimated with 
different theories for the equidistant level case. The exact (red 
solid line), QRPA (black dashed line) and P-QTDA (green 
filled circles) are shown. The ppRPA is also shown (dotted 
line) up to G = 0.48 MeV. 

To further quantify the predictive power of the P- 
q TDA (gs) m ethod, some ratios of pair transfer proba- 
bilities are shown in Figure[7] In this figure, P gs , P ex and 
Ptot correspond respectively to the probability to transfer 
to the ground state, the sum of probabilities to transfer to 
any excited states, while Ptot = Pgs+Pcx- Below the pair- 
ing threshold, the P-QTDA( GS ) reduces to a ppTDA( GS \ 
where the diagonalization is made in a reduced space of 
Slater determinant. In that case, the result are of the 
same quality as for the ppRPA. This figure clearly con- 
firms that while QRPA is rather far from the exact re- 
sults in the superfluid phase, especially in the vicinity of 
the BCS threshold, the projected theory provides a much 
better reproduction of ground state to excited state pair 
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transfer probabilities. Moreover, theoretical predictions 
based on the QRPA [38] might significantly overestimate 
the GP V cross section, as suggested also in j3] , especially 
in the vicinity of the normal-supcrfluid transition. 

V. CONCLUSIONS 

In this work, the QRPA description of the two-particle 
transfer mechanism is tested against the exact solution 
in the Richardson model for several conditions of pair- 
ing interaction strength and level spacing. It is seen that 
both ppRPA in the normal phase and QRPA in the su- 
perfluid region are able to grasp the gross feature of the 
pair transfer process. However some differences are ob- 
served. At variance with other kind of resonance states 
mainly build in terms of particle-hole excitations (see for 
instance Fig. 1 of Ref. [37]]) the particle number con- 
servation seems to play an important role when the par- 
ticle number change during the physical process under 



interest. A method is proposed here to improve the de- 
scription of pair transfer in finite superfluid systems. The 
new method is based on the direct diagonalization of the 
Hamiltonian in a reduced space of the projected ground 
state plus two quasi-particle states. This theory improves 
considerably the description of the pair transfer process. 
On the practical side, the P-QTDA^ GS ' requires only to 
solve the BCS or HFB problem in the initial nucleus and, 
except the additional numerical cost of projection, it does 
not need more effort than the original QRPA. Work is in 
progress to apply it to nuclear transfer reactions. 
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